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The observed flat rotation curves of galaxies require either the presence of dark matter in Newto- 
nian gravitational potentials or a significant modification to the theory of gravity at galactic scales. 
Detecting relativistic Doppler shifts and gravitational effects in the rotation curves offers a tool 
for distinguishing between predictions of gravity theories that modify the inertia of particles and 
those that modify the field equations. These higher-order effects also allow us in principle, to test 
whether dark matter particles obey the equivalence principle. We calculate here the magnitudes of 
the relativistic Doppler and gravitational shifts expected in realistic models of galaxies in a general 
metric theory of gravity. We identify a number of observable quantities that measure independently 
the special- and general-relativistic effects in each galaxy and suggest that both effects might be 
detected in a statistical sense by combining appropriately the rotation curves of a large number of 
galaxies. 

I. INTRODUCTION 

The rotation curves of galaxies are direct probes of the shape of their gravitational fields and their matter content. 
That the inferred circular velocities remain approximately constant, even at large distances from the central luminous 
matter, is the strongest evidence for the presence of dark matter at galactic scales [T]. 

Despite its simplicity and remarkable success in accounting for observations over a wide range of scales, explaining 
galactic rotation curves with dark matter remains a hypothesis; to date the candidate dark matter particles have eluded 
direct detection [5]. In the meantime, several possible modifications to the theory of gravity have been explored in 
attempts to explain the observed galactic rotation curves The MOdified Newtonian Dynamics (MOND) framework 
has been the most successful attempt phenomenologically [3j[4|, but suffers from the fact that it is not relativistic. As 
a result, it cannot be used in its empirical form to generate predictions for gravitational lensing or for the dynamical 
evolution of systems at scales comparable to the Hubble scale. 

Earlier [5] attempts to develop a relativistic theory of gravity that mimics the MOND phenomenology often faced 
fundamental difficulties, such as problems with causality or inadequacies in accounting for gravitational lensing in 
galaxies. More recently, new models have been developed to resolve these difficulties, but they introduce several 
additional fields and auxiliary functions. Such additions negate the most appealing aspect of the original MOND, 
i.e., that galactic rotation curves and the Tully-Fisher relation were accounted for with the introduction of a single 
acceleration scale. 

Most previous attempts aim to reproduce the MOND phenomenology by modifying the general relativistic field 
equations. However, in principle, the MOND phenomenology can also be achieved in the non-relativistic limit by 
modifying the equivalence principle, i.e., the inertia of test particles 7 .Which of the two aspects of the theory of 
gravity need to be altered in order to account for the observed rotation curves of galaxies, in the absence of dark 
matter? This question cannot be resolved solely with observations of the non-relativistic Doppler shifts of tracer 
particles. Ideas involving violations of the equivalence principle have been tested empirically, within the dark matter 
interpretation of galactic rotation curves |H] , but not in the framework of modified gravity 

In this article, we aim to address this question by calculating the second order special- and general-relativistic 
corrections to the Doppler shifts of atomic lines used to infer the rotation curves of galaxies. Similar calculations in 
General Relativity, as well as strategies for using such measurements to map the spacetimes of galaxies have been 
reported in previous studies |9j. Here we focus only on the second-order effects and evaluate them using only the 
symmetries of the spacetime. More importantly, we show explicitly that, if the underlying theory of gravity obeys the 
equivalence principle, then the second-order effects can be determined entirely using knowledge of the non-relativistic 
Doppler shifts and without any assumptions regarding the underlying field equations of the theory. As a result, the 
relation between non-relativistic and relativistic effects can be used as a test of the equivalence principle at galactic 
scales, independent of whether dark matter or modified gravitational field equations are responsible for the flat rotation 
curves. 
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II. PARTICLES AND PHOTONS IN GALACTIC POTENTIALS 



The Circular Orbits of Particles 



We begin by assuming that the gravitational potential of a galaxy exhibits a high degree of axisymmetry. We will 
ultimately calculate the second-order Dopplcr shifts and the gravitational corrections assuming that the spacetime 
of each galaxy is asymptotically flat. This approach is formally appropriate only for a galaxy in an otherwise empty 
Universe; in a subsequent section, we correct for the Cosmological redshift of the galaxy. 

Our assumptions motivate our choice of metric 



ds^ = gtt {r, 9) + 25t0 [r, 9) dt d<j> + g„ (r, 9) dr^ + gee (r, 9) d9^ + 



sin^ I 



(1) 



where gtt, gt<p, drr, and gee are undetermined functions of the coordinate radius r and the polar angle 9. The gt^ 
coefRcient is associated with frame-dragging, which we expect not to be significant for galaxies and will be neglected 
henceforth. For comparison, we recall that in General Relativity, the external spacetime of a spherically symmetric 
object is unique and is given by the Schwarzschild solution, for which gt^ = 0, gee — and. 



9tt = 9r 



1 



2GM\ 
rc^ J 



(2) 



Here G is the gravitational constant and M is the gravitational mass of the object. Because of the assumption of the 
validity of the equivalence principle, matter and photons follow geodesies in the spacetime described by the metric ([I]). 



We describe the motion of a massive particle in terms of its 4- velocity = {u , u'', u 



We choose our coordinate 



system so that the orbit of the particles we study will lie on the equatorial plane, i.e., we will set sin 9—1 and - 



0. 



The requirement 



-1 for the 4- velocity of a massive particle leads to the constraint 



gtt {u'f + grr {u'-f + r^ = 



(3) 



For a particle in a circular _orbit at coordinate radius rg 



we require that the radial component of its 4-velocity is 
zero, which makes equation ([3| a constraint on v} and u** This constraint is also true for the turning points in an 
elliptical orbit. What sets a circular orbit apart is the fact that all points in the trajectory are turning points, i.e., 
that 



du" 
dr 



= 0, 



(4) 



which specifies the u* component of the 4-velocity uniquely 



gtt {ufii- 



1 d\n\gtt\ 

2 d\nr 



= -1 



(5) 



We can use expression (jsj) to solve for the u"* component of the 4-velocity. Whatever our theory of gravity is, we 
expect that the outer regions of the galaxies are in the weak field, so we can find an approximation to the desired 
accuracy by expanding our metric away from the flat solution 



gtt = -1 



2 (2) 

e git 



0{ef 



(6) 



where we have introduced e merely as a dummy parameter that allows us to keep track of the expansion order. 
We denote the leading order correction to the metric element as e^, because we are counting orders in terms of the 
expansion of the velocity, which will be proportional to the square root of g^^\ Using relations — (j6|, we obtain 
that the components of the 4-velocity of a particle in a circular orbit are 





^ 1 


f e 




= 0, 




= 0, 




= e 


-1 



idg^^y 



2 l^" 2dlnr 



0{e) 



2 dlnr 



■0{ef 



(7) 

(8) 
(9) 

(10) 
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We emphasize that the components of the 4- velocity of the particle in circular orbit depend only on the value and the 
local radial derivative of the gtt element of the metric. 

For comparison, in a Schwarzschild metric, the non-zero components of the 4-velocity of a particle in a circular 
orbit are given by 

„' = l + 3^ (II) 

^.U^f. (12) 

r \ r J 



B. The Redshift of Photons 



Our next goal is to calculate the trajectories and energy shifts of photons as they propagate from their origin in 
the galaxy to a distant observer. For a photon with 4-momentum fc'' = {k^ , k"^ , , k'^) , there are conservation laws 
that arise from the two Killing vectors f = (1, 0, 0, 0) and 77^ = (0, 0, 0, 1) of the spacetime, namely the conservation 
of energy 

£p = = -9ttk' (13) 

and of angular momentum 

Ip = g^^k^i^- = r^k'^ . (14) 
In its trajectory, the photon experiences an overall redshift and Doppler shift, which is given by 

1+^^^- 'r^'iTi , (15) 

where the subscripts "e" and "obs" refer to the emitter and the observer, respectively. 

Because of our assumption of asymptotic flatness, at the location of the observer the spacetime is Minkowski and, 
therefore, guiTohs) — > —1- Moreover, because we are considering a static observer, its 4-velocity is u'^y^^ — (1,0,0,0), 



and the denominator of the fraction in equation (15) is equal to — Cp. The 4-velocity of the emitting particle is given 
by relations ([7])-([l0]), which after inserting into equation (15) leads to 



In this expression, we have dropped the subscripts 'e' and 'obs' for the emitter and the observer. It is implicitly 
understood, however, that the redshift z is measured at the location of the observer, whereas all the quantities in the 
right-hand side of the expression are evaluated at the location of the emitter. This expression is identical at this order 
to similar calculations based on different assumptions [9]. 

In order to calculate the quantity l-p/s-p, we first discuss the orientation and geometry of the galaxy, the observer, 
and the photon trajectories in flat spacetime. All corrections due to lensing appear as factors of at least order 0{e^), 
which is beyond the order we are considering here. 



C. The Trajectories of Photons 

We first set a coordinate system (see Fig. [T]) with its origin at the center of the galaxy and oriented in such a way 
that the orbits of the emitting particles lie on the x — y plane. The detector of the observer defines a second plane (the 
image plane) at some great distance D and at an angle t^obs with respect to the direction of the angular momentum 
of the galaxy z. We use the axisymmetry of the galaxy to choose the orientation of the x — y axis so that the center 
of the image plane of the distant observer lies on the y ~ z plane. 

We then set a new coordinate system (indicated by primed quantities) by rotating the original coordinate system 
around the x— axis by angle "dohs- The x' — y' plane of the new coordinate system is parallel to the image plane and 
the unit vector to the image plane is parallel to the z' axis 

z'= (sindohA . (17) 

V COS dohs J 
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FIG. 1: A geometric representation of the galactic disk (grey ) plane P. Photons with momentum k are emitted within a ring 
(dark yellow) about the galactic center and arrive to the observer's image plane P' with impact parameter 6. The image plane 
is at some distance D from the galactic plane, at an angle i^obs with respect to the axis of symmetry. 



Coordinates on the image plane and coordinates on the galaxy plane are related via the rotation 

/I 

= COSlJobs 

V sin ??obs 




(18) 



For an emitter in a circular orbit on the galactic plane at a radius re and at an azimuth with respect to the x— axis, 
this relation becomes 



•1 

cos 1?obs 

, sin i?obs 







cos-i? 



obs 




re COS (pe 

Te COS 'dohs sin ( 
, re sin i?obs sin (, 



(19) 



If we do not consider the bending of the trajectory of a photon due to gravitatonal lensing, then the x' and y' 
coordinates calculated with the last relation will correspond to the location on the image plane where the photon 
emitted by the orbiting object will be detected. The impact parameter of that photon will, therefore, be equal to 



1 - sin^ d„ "="2 



^obs sm 



(20) 



For a photon with wave vector k to intersect the image plane at a right angle and with an impact parameter h we 
require k \\ so fc = \k\ z' = \k\ [smdoy + cosi^o^)- In spherical polar coordinates, for a photon emitted in the plane 
= ^ , the transformation from cartesian to polar is straightforward and gives 

) ■ (21) 



fc = I fc| (^sin -dohs sin (fyer + cos i^obs ^ + sin -dohs cos ( 
The null property of the photon 4-momentum requires 



|fc| =£p + 0(e2) 



(22) 



which leads to 



i0e + O(e') 



(23) 



D. Relativistic Redshifts Due to Galactic Rotation 



We are now in position to calculate the Doppler shift for an axisymmetric potential and examine the properties 
particular to our study of galaxies. From equation (16) and (23) we write 



ezi sin iJobs cos i 



e^Z2 ■ 



(24) 
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where 



Z2 



2 d\nr 



9i? 



2 d\nr 



(25) 
(26) 



Combining the first two orders, we find 



Z2 



Z + 



(2) 

9tt 



(27) 



This last relation expresses simply the fact that the second-order energy shift has two distinct contributions: one 
from the second order special-relativistic Doppler shift (captured by the first term in the above sum) and one from 
the gravitational redshift (captured by the second term in the above sum). 

In a Schwarzschild spacetime, the expressions for the Doppler shift plus redshift to all orders become 



Zl 



Z2 



'GM 
3GM 



1/2 



2r,c2 ' 



(28) 
(29) 



where we introduced appropriate powers of the speed of light c for completeness. Because the Schwarzschild spacetime 
has a single scale, the relation between the first- and second-order terms is the quadratic 



Z2 



(30) 



and second-order effects are always suppressed compared to the first-order effects. 

Observations [TU] suggest that, to leading order, the velocity profiles of galaxies are very nearly flat over the radii 
of interest, i.e., 



Zl 



(^ 



(7)' 



(31) 



where uq is the inferred, nearly constant rotational velocity, at a characteristic length scale tq and we have introduced 
the small parameter a ~ to describe weak deviations from a constant velocity profile. Such a rotation law implies 
that the spacetime of the galaxy is described, to leading order, by 



(2) 

9tt 



In 



K r ) 
I ''out 



2a 



if a > 
if a = 



(32) 



In this last expression for the special case a = 0, we have introduced as an integration constant the radius rout 3> r at 



which the correction term gl^ drops rapidly to zero. If a dark matter halo with a density profile p • 



is responsible 



for the flat rotation curve of a galaxy, then Tout is the outer cut-off of the halo, which is necessary for the total mass 
of the halo to be finite. The exact value of this constant does not affect the calculation of the redshift to flrst order. 



Using equation (32 ) for the it— element of the metric, we now calculate the next order correction to the redshift as 



Z2 = 




if a< 
if a = 



(33) 



For the case of a perfectly flat rotation curve (i.e., a = 0), the sign and magnitude of the second-order wavelength 
shift Z2 depends explicitly on the cut-off radius rout • To avoid the additional complications introduced by the presence 
of the cut-off radius, hereafter, we will assume a 7^ and not discuss any longer the singular case of a perfectly flat 
rotation curve. Surprisingly, the second-order effects can become significant for sufficiently flat rotation curves (i.e., 
when a w 0). Formally speaking, our expansion is valid only for a > wq/c. 

Note that the relation between the first- and second-order effects for a galaxy with a nearly fiat rotation curve is 



1 



Z2 



2a 



1 

2^ 



zl 



(34) 
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Equations ([t]), ( [To| ), (32), and (33 1 represent the main result of the last two sections; that both the velocity profile 
of matter and the gravitational and Doppler redshifts experienced by photons depend only on the same two local 
properties of the metric at the place of emission: the value of its tt— element and its radial derivative. Neither of the 
two quantities depend on the field equations of the theory of gravity, i.e., of the equation that determines the metric 
elements given a distribution of matter. 



III. DOPPLER AND GRAVITATIONAL CORRECTIONS TO THE LINE PROFILES 



We now examine the implications of the results derived in the previous section for the atomic line profiles detected 
from galaxies. The flux an observer detects at a great distance D is proportional to the integral over the image plane 
of the specific intensity of rays that arrive perpendicular to the image plane, i.e.. 



Feie) j dx' j dy' h{e) .x 



(35) 



Using the system of equations ( 19 1 we convert this integral into one over coordinates in the galactic plane 

i)Je{£e,re,<j)c) ( — ) , (36) 



FM = / r^dr^ 



L>2 



where we have used the invariance of the quantity /g/e'^ to relate the intensity at arrival to that at emission by 

3 



(37) 



The observed energy e and the emitted energy Sc are related by the redshift relations derived in the previous section 



= 1 + z(ro, 0o) = 1 + 21 (re) siniJobs cos 0c + 22(^0) , 



(38) 



where we have explicitly denoted the dependence of the redshift experienced by each photon on the location of its 
emission. 

We now assume that the emission at the local Lorentz frame is mono-energetic, at a rest-frame energy Sq. In other 
words, we assume that 



4, (£e, re-, 4>e) = I{rc-, 4>c)5 [ee(''e, '/>e) " ^o] = li^c, '/>o)<5 {s [1 + z{re, (/fe)] - £o} 

The flux integral, therefore, becomes 

cos -dohs 



Feie) 



1)2 



dr^rc d4>eIirc,4>c)S - eo[l + z{re,4>e)] [I + z{rc, 4>c)] ^, 



(39) 



(40) 



where we have made a change of variables in the (5-function to reflect the fact that the right-hand side of this equation 
is a flux density in the observed energy s. We use the (5-function to evaluate the integral over 0o using the relation 



d9_ 
d(j) 



where is each solution to the equation g{(f>i) — or, in our case, 



9( 



e - eo [1 + z{re, (j)i 



Using the expression (24) for the redshift, we obtain 



1 



cos ( 



zi sin dohs ^ £ 



-1-22 



(41) 



(42) 



(43) 



which leads to two solutions for the angle 4>i with opposite signs. Requiring that | cos < 1 allows us to place limits 
on the range of photon energies that contribute to the line as 



1 



1 ± zi sin dohs + Z2 



(44) 
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In evaluating the integral, we also need the derivative 



dz 




d^e 


<t>.='t>. 



Using the above expressions, the flux becomes 



COSt^obs i 
Pe{e) - j 


^drerc^ J(ro,4) [1 + z(rc,( 




^COS l^obs 


j dvcTc [1 + z{rc, (t)t)] ^ 


dz 


-1 




d(l)e 





dg 



1 1/2 



Z2 



(45) 



= 2 



cos do 



I (re, (pi) +I(ro,-(/ii) 



bs 



cos t^obs 



dr. 



dr. 





-1 






d0e 




2 



zl SIT? do 



bs 



£0 - e 



Z2 



-1/2 



.(46) 



To explore the properties of the second-order corrections to the line profiles, we make for now the simplifying 
assumption that the emission comes from a single annulus in the galactic disk with intensity that is independent of 
azimuth, at a radius tq and with a width Sr. In this case, the line profile becomes 



Fs{e) = 2cosi9obsl(fo) 



sin dohs 



£o 



£0 



Z2 



-1/2 



(47) 



which we express in terms of the dimcnsionless quantity 



r{)5r 

2cos'dohJ^{ro) \ L)2 

— z-y sm "dohs - 
£0/ 



£o 



£o 



^2 



-1/2 



(48) 



The minimum and maximum photon energies for which the flux is non-zero are given by equation (44). Because of 
the second-order Doppler shifts and the gravitational redshift, the center of the broadened profile is displaced from 
the rest energy of the line by an amount equal to 



£ = ^ (£_ + £+) = £o (1 + zl SIV? -dohs - Z2) 



(49) 



Moreover, the amplitude of the red wing of the line is smaller than the amplitude of the blue wing. 

Figures ([2]) and (|3| show the dimcnsionless line profiles from narrow annuli in a Schwarzschild spacetime and in the 



spacetime of a galaxy with a nearly flat rotation curve, for different values of the parameters zi and a (see equation 30 



and [34]) and for different inclinations of the observer. In the case of a galaxy with a nearly flat rotation curve, the 
broadened line is gravitationally redshifted by an amount that is larger compared to the equivalent Schwarzschild 
case. If the flat rotation curve is a result of a dark matter halo, the additional redshift occurs because the dark matter 



halo has a density profile p 



and a large amount of mass (and hence of gravitational redshift) exists outside 



the location of the annulus. If this second-order effect is not corrected for, it will be assigned to the overall recession 
velocity of the galaxy Vg that is due to its peculiar motion and to the Hubble flow. The error is, nevertheless, small 
as it is of order 







£ - £0 








£0 





• 2 q 1 +" 

sm iJohs 7. 

la 



Z^ 



1 

2a 



c 



1 fVrr 



ac 



c 



< 



1 /K, 



< 1 



when a > 



(50) 



where T^iot is the inferred rotational velocity from the non-relativistic Doppler shift. 
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FIG. 2: Profiles of broadened atomic lines tliat originate in an annulus in a Schwarzsdiild spacetime for different values of the 
non-relativistic Doppler sliift zi and inclinations of the observer. The overall effect of the relativistic Doppler shift and of the 
gravitational redshift is to make the profile asymmetric and to shift it towards lower energies. 




FIG. 3: Profiles of broadened atomic lines that originate in an annulus in a galactic disk plotted for different values of the 
parameter a, which measures the degree of flatness of the galactic rotation curve. The red line shows, for comparison, the line 
in a Schwarzschild spacetime with the same amount of non-relativistic Doppler shift. The spacetime that corresponds to a flat 
rotation curve leads to a larger overall redshift as well as to less pronounced blue and red wings for the line. Note that the 
parameter zi used for this plot is unphysically large for a galaxy and was chosen here in order to make the effect visible. 



IV. A STATISTICAL MEASURE OF THE RELATIVISTIC DOPPLER SHIFT AND OF THE 

GRAVITATIONAL REDSHIFT 



In the previous section, we calculated the profile of an atomic line that originates in the equatorial plane of a galaxy 
in a general metric theory of gravity. We showed that the non-relativistic Doppler shift, the relativistic (i.e., second- 
order) Doppler shift, and the gravitational redshift experienced by a photon from its origin to a distant observer 
depend on the magnitude and local radial derivative of the tt-element of the metric, gu, at its origin. We derive here 
quantities that enable the measurement of these second-order corrections. 

We start by considering a single galaxy at a redshift Zg, which accounts for both the peculiar velocity of the galaxy 
as well as the Hubble flow. We will assume that the rotational profile of the galaxy, as inferred from non-relativistic 



Doppler shifts, is nearly flat and described by equation (31 ) from some inner radius tq to an outer radius 3> tq, i.e.. 



Zl 

The relativistic Doppler shift and the gravitational redshift introduce an additional overall change in the energy of 



the photon that we described by the quantity Z2 given, in general, by equation (27 1. If photons and particles follow 
geodesies in the same spacetime, i.e., if the theory of gravity obeys the equivalence principle, then the gravitational 
redshift can be calculated using the same metric that determines the velocities of the emitting hydrogen atoms. In 



this case, Z2 is related to Zi according to equation (34 1. If, on the other hand, the theory of gravity does not obey 
the equivalence principle, then the amount of gravitational redshift experienced by each photon will not have the 
same relation to the non-relativistic Doppler shifts. We parametrize the possibility of the theory not obeying the 
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equivalence principle by the single constant / and write, in general, the second-order energy shift as 



Z2 



2a 



(52) 



For a gravity theory that obeys the equivalence principle, / = 1. 

The line profile measured by an observer at infinity will extend between the two energies £_ and £_|_ dominated by 
the largest rotational velocities and, therefore, (see eq. [44]) 



£o 
l + z, 



The width of the line is 



Ae = £+ — £_ 



1 + Zg 



2 ( ^ ) sin i9obs 



and is determined by the non-relativistic Doppler shift. On the other hand, the center of the line is 



1 , 



1 + z„ 



sm -Oobs ^ — 

2a 



(53) 



(54) 



(55) 



and deviates from eo/(l + Zg) by second order effects. These two gross properties of the fine profile depend on four 
parameters of the galaxy and on the gravitational theory: the redshift of the galaxy Zg, its rotational velocity uq, the 
inclination of the observer i?obs) and the ratio f /a. In principle, we can perform an independent measurement of the 
redshift Zg of the galaxy using optical lines from the galactic nucleus. For a sufficiently large sample of galaxies, we 
may also use a statistical argument regarding the distribution of inclinations. Hence we can use precise measurements 
of both Ae and e for a large number of galaxies in order to measure statistically the parameter f /a and, hence, 
constrain deviations from the equivalence principle. 

Because the above argument is of a statistical nature, it will be difficult to determine the formal and, more 
importantly, the systematic uncertainties of the result based solely on this single type of measurement. There is, 
however, an additional measureable quantity that provides an independent measure of the same parameters and can 
be used as a consistency check in case any deviations from the equivalence principle are found. The relativistic- 
Doppler shift and the gravitational redshift not only introduce an additional energy shift to the line, but also make 
it asymmetric, with the blue wing appearing brighter than the red wing. As a result, we obtain an independent third 
observable from each galaxy in the flux averaged photon energy 



Jde e F,{e) 



(56) 



For a symmetric line, (e) = e; any correction to this equality will be due to the second-order relativistic effects. 

For simplicity, we first perform the calculation of (e) in the frame of the galaxy and add the redshift due to the 
Hubble flow and its peculiar velocity only in the final result. In deriving the fiux averaged photon energy, we need to 
evaluate two integrals of the form (see eq. [46]) 



En = Jdee^F,{s) 

cos dohs 



We write 

and make the change of variables from e to 

de = — £o 



dr^rc / de e" [1 -I- z(re, (f)i)] 



dz 


-1 






d(t>e 




2 



£ = £0 



1 



1 + z{rc,(t)i) 



1 



n 2 



1 + z(rc,0i) 



dz 



d(j)i 



(57) 



(58) 



(59) 
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As a result, the general expression for each integral becomes 



Er, 



ccmJoi 



-3—n 



, cos 



obs 



£)2 

,COSJ?obs 
L>2 



d(f)i [1 + z(ro,0j)] 
dr^rc J d(/)jI(ro, [1 + sin?9obs cos (pi + 22] 
rcdrcd(j>il{rc, 0i) < 1 - (3 + ri) 



— 3— n 



(4 + n) 2 2 2 

zi sin i?obs cos 04 + 22 ^ Zi sin i?obs cos 1 



(60) 



Inserting first the general relation (52 1 between the first- and second-order energy shifts, we obtain 



En = £0 



cos 7?o 

~1d^ 



Vedredcpil {re, (pi) + (3 + n) 
f + a 



^ sin^ -dohs I redred(j).il {r^, (pi) cos^ (pi 



2a 



redred(p^T{re,(pi) z^ -sini?obs redred(pil {re,(pi) zi cos (p. 



(61) 



We now write this expression more compactly by defining an intensity average of quantities involving zi and cos (pi 
as follows 



, a B ±\ I redred(pi I (re, (pi) zf COS^ (pi 

{zi COSP (P,) = J n-( —\ 

J redred(pi I (re, (pi) 



The flux averaged photon energy is then given by 



Eo 1 + Zr, 



1 — sin z9obs {zi cos (p) +9 sin^ z^obs {zi cos (p)"^ + 4 sin^ i^obs {zi cos^ 



/ + "/., 2 



2a 



(62) 



(63) 



The fluxed average photon energy in each galaxy depends strongly on the degree of asymmetry in the emission. 
However, we can again make a statistical measurement of the parameter //a using a large sample of galaxies. 
Indeed, because the angle (pi is measured with respect to the distant observer, its values for different galaxies will 
not be correlated with the first-order redshift zi but will be randomly distributed. As a result, in a statistical sense, 



equation ( 63 ) , when averaged over all possible azimuthal orientations of the observer, will be 

2 



El ^ 

Eq l + z„ 



1 



- 



2 sm iJobs 7, 

la 



(64) 



where we have also taken advantage of the slow variation of zi with radius re. When combined with equation (54) 
for the width of the line and averaged over all possible inclinations of the observer 6'obs, the flux averaged energies of 
a large sample of galaxies will also lead to an independent measurement of the parameter //a and hence to a test of 
the equivalence principle. 



V. DISCUSSION 



In this article, we calculated the second-order special relativistic and gravitational effects on the rotationally broad- 



ened line profiles from galaxies. We identified two measurable quantities (eqs. 55 and 64 ), which can be used in 
conjuction with the observed line widths (eq. [54]) to test the validity of the equivalence principle at galactic scales. 

The level at which we can perform this test depends on the formal uncertainties in each measurement, as well as 
on the degree of azimuthal symmetry in the emission of each galaxy, which can mask the asymmtry in the line profile 
due to relativistic effects. To estimate the magnitude of these two sources of uncertainty, we examine some HI profiles 
provided as results of ongoing single beam surveys [llj . 

As evident from Fig. [4] line profiles for individual galaxies, even when they are of high signal-to-noise and sharply 
double peaked, can have morphologies quite different from the azimuthally symmetric case of Fig. [3] This variance 
is because the HI is not necessarily smoothly distributed and galaxies often are lopsided in the stellar or gaseous 
distribution [T^. Therefore, even if the effect we are searching for is large, which it is is not, the nature of this 
measurement must be statistical. 

We also estimate the internal uncertainties using an estimate of the spectral noise from regions outside of the profile 
and 1000 simulations with noise added to the spectra. Because this noise is superposed on the original spectrum, the 
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FIG. 4: HI profiles of two massive galaxies {Vrot > 200 km/sec) that are among those with the highest signal-to-noise ratio 
in the sample. The two profiles show different signs of the asymmetry, reflecting the "noise" introduced in such measurements 
by the internal, structural asymmetries of the galaxies. 

noise of the simulations is actually V2 larger, resulting in a slight overestimate of the internal uncertainties. We find 
that these are always smaller than the systematic errors due to intrinsic profile asymmetries. 

A precision measurement of the parameter /, which measures potential violations of the equivalence principle, 
requires great care in handling of the data as subtle biases can be introduced by careless binning or scaling. While 
we argued earlier that we could average over a suitably large distribution of inclinations, it is, in principle, possible to 
obtain independent measurements of the inclinations of individual galaxies which would further restrict the range of 
the constraints. New surveys will be releasing many thousands of galaxy spectra in the near future, which offers the 
possibility of placing constraints on potential equivalence principle violations on the scale of galaxies. We will discuss 
a detailed observational strategy for dealing with these issues in a forthcoming article. 
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